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ABSTRACT 

The long-standing puzzle surrounding the statistical mechanics of self-gravitating systems 
has not yet been solved successfully. We formulate a systematic theoretical framework of 
entropy-based statistical mechanics for spherically symmetric collisionless self-gravitating 
systems. We use an approach that is very different from that of the conventional statistical 
mechanics of short-range interaction systems. We demonstrate that the equilibrium states of 
self-gravitating systems consist of both mechanical and statistical equilibria, with the former 
characterized by a series of velocity-moment equations and the latter by statistical equilib- 
rium equations, which should be derived from the entropy principle. The velocity-moment 
equations of all orders are derived from the steady-state collisionless Boltzmann equation. 
We point out that the ergodicity is invaUd for the whole self-gravitating systems, but it can 
be re-established locally. Based on the local ergodicity, using Fermi-Dirac-like statistics, with 
the nondegenerate condition and the spatial independence of the local microstates, we red- 
erive the Boltzmann-Gibbs entropy. This is consistent with the validity of the collisionless 
Boltzmann equation, and should be the correct entropy form for collisionless self-gravitating 
systems. Apart from the usual constraints of mass and energy conservation, we demonstrate 
that the series of moment or virialization equations must be included as additional constraints 
on the entropy functional when performing the variational calculus; this is an extension to the 
original prescription by White & Narayan. Any possible velocity distribution can be produced 
by the statistical-mechanical approach that we have developed with the extended Boltzmann- 
Gibbs/White-Narayan statistics. Finally, we discuss the questions of negative specific heat and 
ensemble inequivalence for self-gravitating systems. 

Key words: methods: analytical - galaxies: kinematics and dynamics - cosmology: theory - 
dark matter - large-scale structure of Universe. 



1 INTRODUCTION 

It is known that conventional methods for the statistical mechanics 
of short-range interaction systems cannot be used directly to 
study long-range self-gravitating systems. Hence, it is necessary 
to return to the starting point of statistical mechanics and to 
develop special techniques t o handle the long-range nature of 
gravity jPadmanabbMil Il99(]|) . However, for more than several 
decades, the statistical mechanics of collisionless self-gravitating 
systems has not been realiz ed successfully, and it still remains 
a puzzle (see, for examp le, lOgorodnikov igS^; 'Antonov"l962'; 
Lvnden-Beil Il967l; IShij 119781: Tremairie. Henon & Lvnden-Bell 



1987 



19861 : ISricfcatI Il987l: IStiavelli & BertinI Il987l: Iwhite & Narayan 
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'1992'-" 
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1998 



KuU. Tre umann & Bohringer 1997; Chavanis & Somme ri^ 

Nakamura 2000; Hansen et al.l I200A (Levin. Pakter & Rizzatol 
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In a previous work dHe & Kand l201(t) . we performed a 
preliminary investigation of the statistical mechanics for self- 
gravitating systems. We employed a phen omenological entropy 
form of ideal gas, first proposed by [White & NaravanI (Il987h . to 
revisit this question. By calculating the first-order variation of the 
entropy, subject to the usual mass- and energy-conservation con- 
straints, we obtained an entropy stationary equation. Combined 
with the Jeans equation, and by specifying some functional form for 
the anisotropy parameter /?, we numerically solved the two equa- 
tions. We demonstrated that the velocity anisotropy parameter plays 
an important role in attaining a density profile that is finite in mass, 
energy and spatial extent. If combined again with some empirical 
density profile from simulations, our theoretical predictions of the 
anisotropy parameter, and the radial pseudo-phase-space density 
p/(Jr in the outer regions of the dark matter halo, agree very well 
with the data of A'^-body simulations. The predictions are also ac- 
ceptable in the middle regions of dark halos. Disagreements occur 
in the inner regions, probably because some important physics has 
been ignored. 
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The second-order variational calculus reveals the seemingly 
paradoxical, but actually complementary, consequence that the 
equilibrium state of self-gravitating systems is not only the global 
minimum entropy state for the whole system ~ under the condition 
of energy conservation and the assumption that p is an unknown 
variable function, whereas the two second-order velocity moments 
fyVr and pvf are known invariant functions of the radius r - but 
also, simultaneously, the local maximum entropy state for every 
and any small part of the system. This finding suggests that equi- 
librated self-gravitating systems are not in the maximum entropy 
states, and the local statistical properties of self-gravitating systems 
should be different from their global properties. 

In a follow-up work (Kane & He 201 1), by introducing an ef- 
fective pressure instead of a radial pressure to construct the specific 
entropy, we used the entropy pr inciple. We proceeded in a simi- 
lar way as in iHe & Kand J201(t) and obtained a new entropy sta- 
tionary equation. An equation of state for equilibrated dark haloes 
is derived from this entropy stationary equation, from which the 
dark halo density profiles can be obtained. We also derived the 
anisotropy parameter and pseudo-phase-space density profile. All 
these predictions agree well with numerical simulations in the outer 
regions of dark haloes. This work provides further support to the 
idea that statistical mechanics of self-gravitating systems is feasi- 
ble. We believe that our findings might provide crucial clues for 
the development of statistical mechanics for both self-gravitating 
systems and other long-range interaction systems. 

However, our treatments suffer from the following drawbacks, 
(i) The entropy forms are only chosen phenomenologically, with- 
out robust physical origins. (2) We calculate the first-order variation 
only with respect to the density profile, but we treat other variables, 
such as the pressure, as known functions. Thus, the variational cal- 
culus is incomplete. (3) We could still miss some crucial physics. 
Nevertheless, from these works, we can still draw the following 
conclusions. 

(i) Both the concept of entropy and entropy principle are 
still valid for self-gravitating systems. This finding implies that 
the entropy-based statistical mechanics of self-gravitating sys- 
tems is also valid, and we do not agree that a dynamical the- 
ory can account for the c oarse-grained phase-space density (cf. 
lArad & Lvnden-Bellll2005h . 

(ii) Entropy is extensive. We first constructed the specific en- 
tropy, and then obtained the total entropy in the additive and hence 
extensive way, as in equation (4) of iHe & Kand (201C|) or equa- 
tion (6) of lKang&Hjt201lh . The fact that such an extensively 
defined entropy works very well indicates that a non-extensive en- 
tropy, such as the entropy o fifsaiiis (^1988"), is not necessary. 

(iii) The equilibrium sta tes also consist o f mech anical equilib- 
rium. In lHe& Kanl j20ld) and lKang & h3 (201 ll) . we combined 
the entropy stationary equation with the Jeans equation to predict 
the final state of the system. This proved to be a successful ap- 
proach. The success of this treatment indicates that the final states 
are not only statistical equilibrium states, but also mechanical equi- 
librium states. 

(iv) Local and global equilibria are different. That is, for an 
equilibrated self-gravitating system, the local statistical properties 
differ from the global properties. This difference suggests that we 
should develop the statistical mechanics of self-gravitating systems 
in a different way from the standard statistical mechanics of short- 
range interactions. 

We use the above four principles to formulate the theoreti- 
cal framework of the statistical mechanics for spherically symmet- 



ric coUisionless self-gravitating systems. The paper is organized as 
follows. In Section |2l we present the series of moment equations, 
which are derived from the spherical coUisionless Boltzmann equa- 
tion (CBE) and are employed to characterize the mechanical equi- 
libria of the system. In Section|3] we demonstrate that ergodicity is 
invalid globally for self-gravitating systems, but that it can be re- 
established locally if the particles are treated as indistinguishable. 
In Section|4l we present the entire framework of the statistical me- 
chanics, with a truncated distribution function (DF) to show how it 
is possible to work with this statistical mechanics. We give a dis- 
cussion and conclusions in Section[5] 

2 RELAXATION PROCESSES AND 
MECHANICALA'IRIAL EQUILIBRLV 

2.1 Relaxation Processes 

In physical dynamics, relaxation refers to the essential process 
by which a system approaches equilibriu m or by which a per- 
turbe d system returns to equilibrium (Mo, v an den Bosch & White! 
I2OIOI) . So far, four coUisionless relaxation processes have been 
identified as being responsible for the equilibrium state of self- 
gravitating systems: (i) violent relaxation, (ii) Landau damping, 
(iii) phas e mixing and (iv) chaotic mixing. We refer the reader to 
lMoetalJ ( l201Q) for a detailed discussion of these relaxation pro- 
cess es. There could also be other possible relaxation mechanisms 
(e.g. iGurzadvan & Savvidvl[T986l) . 



2.2 Velocity moment equations and mechanical equilibria 

The kinetic-theoretical description of self-gravitating systems is 
well approximated by the CBE as 

1^ = ^+v^- — -^=0 (1) 
At dt dx dx dv 

Here, f{t,x,v) is the fine-grained statistical DF. The CBE is an 
approximation that results from the Liouville equation when the 
particle number of the system is very large, and the A^'-body DF 
can be separated as the product of single-particle DFs as 

N 

{t,Wl,...,WN) = Y[f{t,W,). (2) 

Here, Wi = {xi,Vi) denotes the phase-space coordinates for the 
i-th particle. This assumption implies that positions of gravitating 
particles are uncorrelated (i.e. the probability of finding a particle 
near any phase-space position is unaffected by the pre sence or ab - 
sence of particles at nearby points; see iBinnev & Tremainell2008h . 
We refer to this separability of the A^-body DF as the statistical 
independence of particle positions. Furthermore, the gravitational 
effect exerted on a specific particle is collectively represented by 
the mean gravitational potential ^{x), contributed from all the par- 
ticles of the system. This mean background potential <I>(jc), as well 
as the uncortelation assumption, suggests that the direct gravita- 
tional interaction between any two particles is unimportant, and 
hence should be neglected (i.e. the particles among themselves are 
coUisionless). 

The CBE might not be a good approximation in very dense re- 
gions, such as the center of a dark matter halo, where two-body re- 
laxation might play an important role. In this case, the CBE would 
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be be t ter replaced by, say, the Balescu-Lenard equatiorQ ( lHevvaert3 
bOldlChavaiiislbOlOt) . 

Throughout this work, spherical symmetry is always assumed 
for self-gravitating systems. However, we believe that it is not diffi- 
cult to extend our results to general cases. We multiply the steady- 
state CBE with various powers of velocity, and we integrate over 
all velocities to obtain the moment equations of all orders (see Ap- 
pendix [A^. The second-order equation is the familiar Jeans equa- 
tion: 



d , 2p,^ d , 28 d$ 

dr r dr r dr 



(3) 



Here, Vt denotes either vg or and the barred quantities indi- 
cate the corresponding velocity moments. /? is the usual velocity 
anisotropy parameter, defined sls 13 = 1 ~ We can also ob- 

tain two 4th-order equations as 



d_ 
dr 
and 
d 



{pvt) + ^{vt- 



4p 



- 3v?.v? 



'3pVr 



rd$ 

dr ' 



(pv?v?) + ^{3v?v^ ^ vf) 



These two equa t ions are actually the same as th ose of 

iMerrifield & Kenll jl99(]|) . except that IMerrifield & KentI jl990l) 
used Vt to denote the two-dimensional velocity component, as 
Vt = {vg + v'^Y^^ . The general moment equation s are give n by 
equations iATB (see also Deionghe & Merritt.l99a : iAr]|201lh . 



Td$ 



(4) 



(5) 



dr 



{pVr 



(m + 2) 1, 1 5 — 
+ — —^pv'^+^v]; 



{k + 1) (m + 2) 

r (m + 1) 



X pv^vr^^ 



-{k- 



,d£ 
dr ' 



l)-:—pvpvY^. 



Generally, there are A*' independent 2A'^th-order moment equations, 
with different combinations of k and m, in equation dAl It . These 
moment equations characterize the complex mechanical equilibria 
of self-gravitating systems. 

As demonstrated in Appendix [A] the DF / is insepara- 
ble, and hence these infinitely many moments are independent 
and irreducible. The right-hand side of equation (lAl U . — (fc + 
l)pi'rti"M$/dr, indicates the self-coupling between the system 
and its own gravitational potential through various orders of veloc- 
ity moments specified by different k and m. 

Magorrian & Binnev ( 1994) have shown that when moments 
up to the lOth-order have been determined, accurate predictions of 
the line-of-sight velocity distribution can be made. This result sug- 
gests that a complete and accurate prediction of the 3-dimensional 
DF might also need to be calculated to this high order. 

In Fig. [T] we show independent velocity moments up to 
the lOth-order for spherically symmetric self-gravitating systems. 
These non-vanishing moments are the effective thermodynamic 
variables for self-gravitating systems. Obviously, the number of 
equations of mechanical equilibria (i.e. the moment equations lAl lb 
is less than the number of independent moments. The series of mo- 
ment equations dAl It does not constitute a closed equation sys- 
tem, because the moments pv^, the one-dimensional transverse 
moments of all orders (fc = 0,2,4,...), do not appear in the differ- 
ential terms, and hence there are no differential equations for them. 

The DF / that satisfies the CBE is the fine-grained DF, or 
phase-space density. Any physical quantity constructed with the 



How ever, the two-body relaxation was overemphasized in iHe & Kand 

iioiol) . 



Order 



■4vf 
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pVr 


pvl° 


pvfvl 


pVrVf 


pVrVf 


pvfvf 



pvf ^-^^VrVf pvfvf 



Figure 1. Independent velocity moments up to the lOth-order for spheri- 
cally symmetric self-gravitating systems. The moments in the upper-right 
part are govemed by the moment differential equations l lAllK i.e. the mo- 
ments appear in the differential tenns. However, the moments pv^ in the 
diagonal, with fc = 0, 2, 4, . . ., do not appear in the differential terms, and 
hence there are no differential equations for these. 



fine-grained DF /, s uch as the entrop y S{f), should be conserved 
(i.e. d5/dt = 0; see lMo et al.ll2010h . The meaningful DF for sta- 
tistical mechanics i s the so-called coar se-grained DF fc(x, v), usu- 
ally d efined as (see lSaslawl2"000l : lBindoni & SeccQ.200a : ,Mo et all 
l20ich 



(6) 



Mx,v)^-^[ /(x',v')dVdV. 

This is the average of / within a sufficiently small but not infinites- 
imal phase-space element, i.e. the macrocell, A^vf — A'^jcA'^v, 
centered at the phase point (x, v). By definition, fc always satisfies 
fc ^ /. Also, it is not difficult to prove that the statistical indepen- 
dence suggested by equation (|2} is also held for the coarse-grained 
DF 

With the definition of the coarse-grained DF, we can find the 
relationship between the coarse-grained velocity moment Qc and 
the fine-grained moment Qf of any quantity Q{v), a smooth func- 
tion of V. We have 



1 



A**vf 



0(v')/(x',v')dVdV~Q(v)/e(x,v), 



(7) 



where Q(v') ~ Q{v) within the small phase-space element A^w, 
and we use the definition of the coarse-grained DF of equation 
Integrating both sides of equation (|7} over v, we obtain 

Q:ix) = j Q{y)fc{x,v)d\=^ j^^ Q(v')/(x',v')dVdV 



/S?x 



A^X 



Q-f{x')d\' ^ Qjix) + ^ ^ ^^Ax?, (8) 
1=1 



where the last approximation is obtained by Taylor expansion with 
respect to x to the second order. Hence, we can see that, although 
the coarse-grained DF does not obey the CBE, the coarse-grained 
moments differ from the fine-grained moments by only up to the 
second-order small quantities determined by the size of the finite 
volume element. So, we can safely use these moment equations to 
characterize the mechanical equilibria of the system, regardless of 
whether these moments are derived from the fine-grained or coarse- 
grained DF. 
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Figure 2. The virialization ratio vs. the radius r of a virialized cold dark 
matt er halo, which is calcu late d using equation lll2t . The simulation data are 
from lNavarro et allteOlOl) and lLudlow et al.lfcOlOl) . The result also implies 
that the CBE is a very good approximation for describing self-gravitating 
systems, such as dark matter halos. 



2.3 Generalized virial equations and virial equilibria 

Multiplying both sides of the velocity-moment equations dAl It by 
A-Kr '"^ and integrating over the radius r, we can obtain the general- 
ized virialization relations as 



-2(m^ - l)rfc+2,™ + 2{k + l)(m + 2)Tk,m+2 
+ {k+ l){rn + l)Wk,,^ = 47r(m + 1^ Pk+2,„. 



in which 

Wk,m : 



(9) 



-d$ 




pVrvV^ ——Aur dr 
dr 



(10) 



Figure 3. Illustration of the ergodicity validity for short-range interaction 
systems and the global ergodicity breaking for self-gravitating systems. 
Panel (a) shows the microstates of a short-range interaction system (e.g. 
an ideal gas), in which the motion states of two particles, as a result of 
stochastic motions driven by irregular forces, can be exchanged with each 
other Panel (b) shows the typical non-resonant regular orbit of a gravitat- 
ing particle in the static spherical potential of the self-gravitating system, in 
which the global ergodicity is invalid. 



andPfc 



pvrvp. When fc = m = 0, equation ^ reduces to 



2fctot(r) + u(r) = 4Tvr Pr(r), 



(11) 



where fctot(r) = r2,o + 2ro,2 andu(r) — Wo,o are the total kinetic 
and potential energy contained in the r-sphere, and Pr(r) — P2,o 
is the radial pressure. By differentiating equation (|9) with respect 
to r, we can immediately restore the moment equations JAllt . It 
can be seen that the moment equations are equivalent to the gener- 
alized virial equations (|9), and hence the mechanical equilibria are 
identical to the virial equilibria. 

From equation Jilt , we define the virialization ratio as a 
function of r as 



rv{r) = 



2fctot(r) 



(12) 



With the density profile p(r), the velocity dispersion profiles cr^fr) 
and o " t (r) of the simula ted dark matter halo (see iNavarro et al.l 
I2OIOI : iLudlow et a we calculate rv, as shown in Fig. |2] It 

can be seen that nearly equals one at all the radii, except for the 
inner fluctuation and a peculiarity at r ~ lOOkpc/h. This agree- 
ment indicates that the simulated dark halo is well virialized, and 
it also implies that the CBE is a very good approximation for de- 
scribing self-gravitating systems, such as dark matter haloes. 



3 ERGODICITY GLOBAL BREAKING AND LOCAL 
VALIDITY 

Ergodicity, or equiprobability principle, is the fundamental for 
equilibrium statistical mechanics. In physics and thermodynamics, 
the ergodic hypothesis says that, over long periods of time, the time 
spent by a particle in some region of the phase space of microstates 
with the same energy is proportional to the volume of this region 
(i.e. that all accessible microstates are equiprobable over a long 
period of timefl Below, we illustrate that the ergodic hypothesis 
is not valid for the whole self-gravitating system, but can be re- 
established in any local finitely small volume element. 



3.1 Complete relaxation and ergodicity validity 

In the usual thermodynamic system with short-range interactions, 
such as an ideal gas, when ergodicity is valid, the constituent parti- 
cles are driven by irregular forces and perform Brownian motions. 
As a result, the particles can traverse any place of the system, and 
can attain any possible value of the momentum (or velocity). In 



^ This statement of ergodic hypothesis is taken from Wikipedia at 
http://en.wikipedia.org/wiki/Ergodic_hypothesis 
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Figure 4. Illustration of the re-establishment of the local ergodicity for collisionless equilibrated self-gravitating systems. For simplicity, the motion states at 
only two instants are shown, (a) Within a finitely small volume element of the system, depicted by the fictitious dashed square, all the collisionless gravitating 
particles have detemiinistic trajectories and labels, and are moving extremely rapidly in randomly oriented straight lines. At the next instant, the original 
particles move out of the volume element, and are superseded by some other particles from the outside, (b) If, within this very volume element, these randomly 
and homogeneously distributed particles are treated as indistinguishable, they will exhibit a static distribution pattern. These patterns at different instants 
constitute an ensemble, every member of which occurs with equiprobability in the hypersuriace, defined by all the invariant velocity moments, for the volume 
element. In this sense, the microstates of the element (i.e. the local microstates) can be defined, and the local ergodicity can be re-established. 



this sense, the stochastic motion of particles and the complete re- 
laxation of the system are, respectively, the necessary and suffi- 
cient condition for the validity of ergodicity. For instance, in an 
ideal gas, the motion states at some instant ti of the two parti- 
cles A and B in Fig.lSja) are and {q2,P2)' respectively. 
Because of the randomness of the motion, at some later time t2, 
the microstate of the system is also available by exchanging the 
motion states between A and B. If the two microstates before 
and after the exchanging of the two particles are equiprobable, i.e. 
-P[^(?i. Pi). -8(92,^2)] = ■P[^(92./'2). -6(91 then the er- 
godicity is valid, and vice versa. 

3.2 Incomplete relaxation and ergodicity global breaking 

However, the above statement of ergodicity is invalid for the whole 
collisionless self-gravitating system. Fig. Ob) shows the typical 
regular, or non-resonant, orbit of a gravitating particle moving in 
an equilibrated spherical gravitational field. An unclosed curve re- 
sembles a rosette in an orbital plane, and eventually passes close 
to every point in the annulus of an invariant torus, on which 
the energy and angular mo mentum of the particle are conserved 
teirmev & TremaiMl2008ll . Even if the orbit suffers from some 
smsill perturbations, according to the Kolmogorov - Arnold - Moser 
theorem, the torus corresponding to the regular orbit with funda- 
mental frequencies sufficiently incommensurable will survive the 
small perturbations. It retains its topology such that the motion in 
its vicinity re mains regular an d confined to the slightly deformed 
invariant torus jMo et al]|2OI0l) . Hence, such regular motion of the 



particle is definitely not ergodic motion, in which the particle can- 
not traverse any place o f the syste m, and also cannot attain any 
value of the velocity (see lHadl2004l) . 

Additionally, Ivan Albadj ( [T982h indicates that violent relax- 
ation cannot wipe out all the information concerning the initial 
conditions when the system settles down into an equilibrium con- 
figuration (see their fig. 3). This correlation between the initial and 
final states also implies that violent relaxation cannot give rise to 
the global ergodicity for the whole system and hence the relaxation 
is incomplete. 

iLvnden-Beiil dlQeTh derived a DF of self-gravitating systems, 
on the basis of the complete relaxation and global validity of ergod- 
icity. However, the Lynden-Bell distribution leads to the isothermal 
mass profile o f a gravitating system, with infin ite mass, energy and 
spatial extent. 'Lynden-Belll l ll967h . lshi] il9J§) and' MadsenI Il987h 
have already suggested that some kind of incomplete relaxation 
me chanism shou ld be considered. For example, the suggested cure 
by iLvnden-Belll i ll 967.) for incomplete relaxation is to introduce 
a cut-off in phase-space, so that complete violent relaxation takes 
place in a limited region only. Shu's approach, however, is to con- 
sider additional macroscopic constraints (Shu 1969, 1987). 

Self-gravitating systems are not the onl y long-range interac- 
tion systems in which ergodicity is broken. ICampa et al.l ( |2009|) 
have reported that many ergodicity-breaking events can be found 
in other long-range systems, such as the generalized Hamiltonian 
mean field mo del, the Ising and th e XY long- plus short-range 
model (see also lBouchet et alj|2008l) . So, it seems that the ergodic- 
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Figure 5. Illustration of the local ergodicity of indistinguishable particles for collisionless equilibrated self-gravitating systems. The upper three scatter plots are 
2-dimensional slices of spatial distributions of the indistinguishable particles within a finitely small volume element. The lower plots are the corresponding 1- 
dimensional velocity distributions within this very element. The local ergodicity is valid, provided that the gravitating particles are treated as indistinguishable, 
that these particles are randomly and homogenously distributed within the finite volume element under consideration, and that the velocity distribution is 
independent of the spatial coordinates. The arrows indicate the time evolution of these distributions for this volume element. We see that, if the possible 
fluctuations are neglected, both the spatial and velocity distribution remain unchanged with time when the local ergodicity is established and the system is in 
equilibrium states. 



ity breaking could be a common feature of long-range interaction 
systems. 

3.3 Incomplete relaxation and ergodicity local validity 

Ergodicity is the basis of equilibrium statistical mechanics for any 
thermodynamic system. If ergodicity is completely invalid, then 
the equilibrium statistical mechanics cannot exist. Fortunately, we 
can re-establish the local ergodicity for equilibrated self-gravitating 
systems in the following way. 

Imagine that, at some instant, within a finitely small (not in- 
finitesimal) volume element of the system, all the collisionless 
gravitating particles have deterministic trajectories and labels, and 
are moving extremely rapidly in randomly oriented straight lines, 
as shown in Fig.|4{a). At the next instant, some of the original par- 
ticles move out of the volume element, but some others will move 
into this very element. As time passes, this process occurs repeat- 
edly and continuously, so that the particles are constantly super- 
seded. Because these particles are not confined within this volume 
element, if the observer's field of vision is restricted within this el- 
ement, these particles cannot trackable. In this case, it is not possi- 
ble to define the microstates of the local element. However, if these 
particles are treated as indistinguishable, even if they have definite 
trajectories and labels, then the physical picture of the microscopic 
motions is renewed, as shown in Fig.lUb). In this picture, at every 
instant, these particles, which are losing their trajectories and la- 
bels, are randomly and homogeneously populated in the volume el- 
ement under consideration. The homogeneous distribution of these 
indistinguishable particles at an instant exhibits a static distribution 
pattern, and these different patterns at different instants naturally 
constitute an ensemble. When the system is in its equilibrium state, 
all the background field quantities (i.e. the velocity moments given 
by equations lA4l orl25t fluctuate at a magnitude proportional to the 
inverse of ^fKN, where AA^ is the particle number of this volume 



element, with /S.N ^ 1. So, these moments can be considered as 
local invariants, which define a confined hypersurface for this vol- 
ume element. It is reasonable to assume that every member of the 
ensemble of the local distribution patterns occurs with equal prob- 
ability in this hypersurface. So, in this sense, the microstates of 
the volume element (i.e., the local microstates) can be defined, and 
hence the local ergodicity can be re-established. 

We provide an intuitive understanding of how to identify the 
local ergodicity. If the particles are treated as indistinguishable, if 
these particles are randomly and homogenously distributed within 
the finite volume element under consideration and if the velocity 
distribution is independent of the spatial coordinates, then the local 
ergodicity can be re-established (see Fig. |5)- Such a requirement 
is extremely easy to satisfy, and hence the local ergodicity should 
always be valid. 

In the previous subsection, we proved the invalidity of global 
ergodicity by regarding the gravitating particles as distinguishable, 
with deterministic trajectories and labels. Even if the particles are 
globally treated as indistinguishable, the validity of global ergod- 
icity cannot be recovered. According to the ensemble theory, the 
probability distribution of the microstates of a closed system at a 
fixed energy is derived by assuming that all states on the energy 
hypersurface in phase spa ce have equal probability. However, as 
Isiimev & Tremain3 (|2008') have argued, the hypersurface of a con- 
stant energy for an isolated self-gravitating system is unbounded, 
and thus the microcanonical probability distribution cannot be de- 
fined. This implies that ergodicity, or equiprobability, is broken 
globally. 



3.4 Negative specific heat and ensemble inequivalence 

The paradox of ens emble inequivalence comes from the fol- 
lowing argument (see iLvnden-Bell & Wood|[l968l : lThirrind[T970l; 
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Figure 6. (a) Illustration of the Fermi-Dirac-like statistics for self-gravitating systems. The local finitely small phase-space element (i.e. the macrocell) is 
divided into many microcells (nine in the figure). The possible minimum phase-space volume is AxAv = hg, whose value is small enough so that each 
microcell can accommodate only one gravitating particle (as it is impossible for two gravitating particles to be close to each other but simultaneously to be 
relatively static). However, the value is large enough that there ai'e no two-body encounters between particles in adjacent microcells. Hence, the microstates 
of the indist ingui shable gravitating particles within this macrocell should be counted by F ermi-Dirac-like statistics. Our approach is more like the particle 
approach of Shu (1978), but such an exclusion effect was already noticed bv lLvnden-Bei] Jl967l) . when he designed his 'fourth statistics'. We assume that 
there are two (indistinguishable) particles. Because of the local ergodicity, the two particles will populate these microcells with equal probability, which also 
implies that all the microstates for this macrocell occur with equal probability. According to the Fermi-Dirac-like statistics, there are in total 36 microstates for 
the two indistinguishable particles in the macrocell (only four of them are shown), (b) Illustration of the spatial independence of the microstates. Here, A and 
B are two macrocells that are spatially adjacent. Because both the direct gravitational interaction and the statistical correlation between any two gravitating 
particles can be neglected, which is exactly the same reason to justify the CBE, it is reasonable to assert that the microstates of the two particles in B do not 
depend on the motion states of the three particles in A. That is, the distribution patterns of the two particles at the two instants in B occur with equiprobability, 
and hence the spatial independence of the local microstates can be established. 



iHeitel & Thirrindll97ll: ICampa et al1l2009l) . In the canonical en- 
semble, the heat capacity at constant volume is 

d{E) 



Cv = 



dT 



= P\{E - {E)Y) > Q, 



(13) 



where {E) is the mean value of the energy. Whereas for the self- 
gravitating system at constant energy (i.e. in the microcanonical 
ensemble), according to the virial theorem, ([/) = ~2{K), so 
E = (K) + (U) — —{K). Because the kinetic energy defines 
the temperature, as T oc A", we obtain the heat capacity of the 
microcanonical ensemble as 



Cv 



dE dE 



< 0. 



(14) 



With the positive and negative specific heat derived in equa- 
tions il3\ and il4l . respectively, we can see that the microcanoni- 
cal ensemble and the canonical ensemble are inequivalent for self- 
gravitating systems. 

We propose that such a paradox should originate from the mis- 
use of the ensemble theory. In fact, a detailed analysis indicates that 
the ensemble theory is established only upon the validity of global 
ergodicity, and hence is only applicable for short-range interaction 
systems. In Section [T2l we have demonstrated that the global er- 
godicity is not valid for self-gravitating systems, so the relation- 
ships of equations J13l > and l ll4t . and any other conclusions based 
on ensemble theories , are rootless for self-gravitating systems (e.g. 
IChavanisll2005l , l2006l) . 

As mentioned in Section [33llBinnev 8i Tremaind ( |2008[) have 
presented an alternative argument against the existence of the mi- 



crocanonical ensemble for self-gravitating systems. They have sug- 
gested that the hypersurface of a fixed energy for an isolated self- 
gravitating system is unbounded, and thus the microcanonical prob- 
ability distribution cannot be defined. Moreover, they have also 
demonstrated that the canonical ensemble does not exist. Hence, 
without the validity of ensemble theory, the relationships of equa- 
tions il3\ and il4\ should be, at most, considered as two different 
definitions of the so-called 'heat capacity' for self-gravitating sys- 
tems. 

In the next section, we demonstrate the significance of this 
local ergodicity in deriving the Boltzmann-Gibbs entropy for self- 
gravitating systems. 



4 BOLTZMANN-GIBBSAVHITE-NARAYAN STATISTICS 

As addressed in the introduction , in the previous studies 
l iHe & Kansl[201cl : lKang & Hdl201ll ). we found the following, (i) 
Both the concept of entropy and the entropy principle are still valid 
for self-gravitating systems, (ii) Entropy is additive and hence ex- 
tensive, (iii) The equilibrium states consist of mechanical equilib- 
ria, (iv) The local statistical properties differ from the global prop- 
erties. Indeed, we have demonstrated in Section [T2l that the global 
ergodicity for self-gravitating systems is broken. However, if the 
particles are treated as indistinguishable, then the local ergodic- 
ity can still be re-established. Such a characteristic of ergodicity 
suggests that the statistical mechanics of self-gravitating systems 
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should be quite different from the standard statistical mechanics of 
short-range interaction systems. 



4.1 Boltzmann-Gibbs entropy and the coarse-grained DF 

We first derive the appropriate entropy form for self-gravitating sys- 
tems, on the basis of the local ergodicity and the statistical indepen- 
dence of particle positions suggested by equation (|2j. 

The local ergodicity implies that the indistinguishable grav- 
itating particles should obey Fermi-Dirac-like statistics (see 
Fig.|6^). We assume that the center of a finite volume element is lo- 
cated at X = {xi ,X2,X'j,), whose volume is A'^jc. The particle num- 
ber in this element is /^N{x) = pjv(x) A^x, where pN is the parti- 
cle number density, a function ofx. By specifying the particle mass. 
Trip, we can convert pjv to mass density pm as pm = m^pN- Let 
n{x,v) denote the velocity distribution, which is the particle num- 
ber in the finitely small phase-space element (i.e. the macrocell), 
A^jcA^v, located at the phase point {x^v). According to the defi- 
nition of the coarse-grained DF fc{x,v) in equation {6}, we have 
n{x,v) oc /c(x, v)A'^xA'^v, and hence AA^(;i:) — y^^, n{x, v). As 
analyzed above, the distribution n{x, v) should obey Fermi-Dirac- 
like statistics. So, the number of the microstates in this volume el- 
ement is 



n 



j{x,v)\ 



n{x, v)l{bj{x,v) — n{x,v))\ 



n 



ojix,v) 



n(x,v) 



n{x, v)! 



AiV! 



AN\ 1 lln{x,v) 



sn(jr,v) 



wmb [n,x\ 
AN\ 



(15) 



where uj{x, v) = A'^xA'^v/hg denotes the degeneracy at the veloc- 
ity level V (not energy level) for the macrocell. If, for all velocities 
V, n{x,v) /ij{x,v) <^ 1 (i.e. if the distribution satisfies the non- 
degenerate condition), we obtain the approximated equality in the 
first line of equation ( I15t . Here, WMBlrijX] denotes the microstate 
number counted with the Maxwell-Boltzmann statistics, given by 
the term enclosed in the parentheses. 

Because the close encounters between gravitating particles 
can be neglected, which is the reason for justifying the CBE, in this 
case hg can be treated as an infinitesimal quantity (i.e. hg — )■ 0), 
and hence ui — ^ oo. As a result, the validity of the CBE is consis- 
tent with the non-degenerate condition, n/cj <^ 1. Whereas in very 
dense regions, such as the center of the dark matter halo, where 
two-body relaxation might play an important role, the CBE might 
not be a good approximation. In this case, the non-degenerate con- 
dition will be invalid. 

As illustrated in Fig. |6jb), it is reasonable to believe that the 
microstates of different macrocells are spatially independent of 
each other. Hence, the total number of the microstates in the whole 
system is constructed as the product of wlujx] of all the volume 
elements. 



W[n] — JJ^ 



(16) 



where TO[n,x] and W[n] are both functionals of the distribution 
n{x,v). Taking the logarithm of both sides of equation ( 116) , we 
have 

In W[n] = — y ^ lnn(j:, v)! + ^ ^ n{x, v) lncj(x, v) 



x,v 



n{x,v)lnn{x,v) + y n{x,v)\nuj{x,v) 



x,v 



(17) 



Here, we use Stirling's formula. Inn! « n(lnn — 1), for the ap- 
proximation in the second line, and A^t is the total particle number 
of the system: 

A^T = ^AN{x) = ^n{x,v). 
X x,v 

By dropping the unimportant constant A^'t, with 



F{x,v) 



n{x,v) 



hi 



A^xA^v 



n{x,v) oc fc{x,v), 



(18) 



we obtain 

In W[F] — — n{x, v) \nn{x, v) + n{x, v) lnu!{x, v) 



E 

x.v 

^ n{x, v) ^ n{x, v) 



uj{x,v) uj{x,v) 



/ F{x,v)lnF(x,v)d^xd^v. 



(19) 



In the last line, ^ A^j:A^v is replaced by J d'^xd^v. Ignoring the 
unimportant factor 1/hg, we use the Boltzmann relation to define 
the entropy, 5 = In W, and then the total entropy of the system is 
exactly the Boltzmann-Gibbs entropy: 



S[F] = F{x,v)lnF{x,v)d^xd'^v. 



(20) 



As mentioned previously, the entropy for self-gravitating systems 
should be additive and hence extensive. We can see that the ex- 
tensivity of entropy comes from the spatial independence of mi- 
crostates (equation [T6l . As in standard statistical mechanics, the 
factor 1/AA'^! in equation dlSt ensures the extensivity of the en- 
tropy of equation ( I17t or equation ( 120) . 

From equation ( 118) , we see that the DF F{x, v) is the same as 
fc{x, v) except for a constant coefficient. Thus, we can regard them 
as identical. The DF F{x, v) can be expressed by Taylor expansion 
with respect to v, as 



lnF{x,v) = 



So, we have 



d 



fc+m + n 



\nF{x,v) I 



^-^ k\m\n\ dv^dv^dvo 

k,m,n 



k m n \ 

V1V2 1^3 -(21) 



F{x,v) 



with 



cxp 



E\ / \ k m n 

y^k,m,n(X)ViV2 V3 



(22) 



d 



lnF(jc,v) 



(23) 



Xk,m,n{x) k\m\n\ dv^dv^dv^ 

where \k,m,n{x) are functions of the spatial coordinates x. We 
should emphasis that the fine-grained DF f{x,v) obeys the Jeans 
theorem that, as a steady-state solution of the CBE, f{x,v) depends 
on the phase-space coordinates only through integrals of motion in 
the given potential, such as the Hamiltonian H, or angular momen- 
tum L. However, the coarse-grained DF F{x, v) does not obey the 
CBE, and hence it is not subject to a constraint of the Jeans theo- 
rem. 
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The formulation of F{x, v), which we have adopted in equa- 
tion i22l . is the fundamental expression. However, the formulation 
is not unique, and we can choose various alternative forms to ex- 
press F{x,v). With the fundamental expression in equation i22i . 
we can define a quantity, Zg{x), as 



Z,ix)= / F-o{x,v)d\ 



(24) 



where F^oix,v) is defined by removing the zero-order term 
from F{x,v), as F-.o{x, v) = F{x,v)/ exp(— Ao,o,o)- The non- 
vanishing velocity moments of various orders, Mk,m,n{x), relate 
to the coarse-grained DF as 



Mk,m,n{x) = P^l^T^Si^) = / vtvTv^F{x,v)d-'v 



exp(-Ao.o.o) / v^vTv^F.o{x,v)d^v 



. ^ . dZg(x) .dlnZg{x) 
= -exp(-Ao.o,o) = -pix) 



(25) 



dA dAfe 

Here, k, m, n are non-negative integers (i.e. k,m,n — 0, 1,2,.. .), 
p{x) is the mass or number density, and we have used 

f , ^ Mofifi{x) p{x) 
exp(-Ao,o,o) = = 

As explained in Section 12.21 the coarse-grained moments differ 
from the fine-grained moments by only up to second-order small 
quantities, determined by the size of the finite volume element. So, 
the moments evaluated by equation l l25t can be directly substituted 
into the moment equations dAl lb . 

We can see that the function Zg (x) greatly resembles the par- 
tition function of the usual statistical mechanics of short-range sys- 
tems. Here, we also call Zg{x) the partition function - of the self- 
gravitating systems. 

From equation ( I25t . we can see that any moment can be de- 
rived from the partition function, Zg(x), by differentiating with re- 
spect to the corresponding coefficient of Taylor expansion. How- 
ever, it would be a great mathematical challenge to evaluate simul- 
taneously both the partition function Zg{x) and the entropy ^[-F] 
of equation l |20t . We leave these mathematical problems for future 
work. 

To summarize, we identify the local ergodicity for self- 
gravitating systems, and we rederive the Boltzmann-Gibbs entropy 
by using Fermi-Dirac-like statistics, the non-degenerate condition 
and the spatial independence of local microstates. The Boltzmann- 
Gibbs entropy is confirmed to be suitable for the statistical mechan- 
ics of self-gravitating systems. 

We emphasize once more that Binney's argument cannot in- 
validate the Boltzmann-Gibbs entropy for self-gravitating systems, 
because we have already revealed that equilibrated self-gravitating 
system s are not in maximimi entropy states (sec iHe & Kana.201Q . 
l201lh . 

4.2 Macroscopic constraints and extended White-Narayan 
statistics 

As mentioned in the introduction, the final relaxed state of self- 
gravitating systems is also the statistical equilibrium state. There- 
fore, it is the task for the statistical mechanics to completely deter- 
mine these expansion coefficients in equation i22l . 

It is easy to think that mass, energy, momentum and angular 
momentum conservation are possible macroscopic constraints for 



constructing the statistical mechanics of self-gravitating systems. 
If the DF is an even function of velocity v and the system is spheri- 
cally symmetric, then the momentum and angular momentum con- 
servation are automatically satisfied, in that they are both vanish- 
ing. Here, we do not consider the case of non-vanishing angular 
momentum. Besides these usual conservation constraints, however, 
the conservation of the phase-spa ce volume of the fin e-grained dis- 
tribution should also be included ( lLvnden-Bellll967h . Moreover, as 
suggested by Shuj (19691, Il987l) . the incomplete relaxation should 
be properly han dled by introducing some macroscopic constraints. 
iMadsed fT987h has also suggested that the Vlasov equation should 
be properly considered. 

We argue that, besides the usual constraints of the given mass 
and energy, the series of moment equations JAl 1) of the steady- 
state CBE must be included as additional constraints on the sys- 
tem's entropy fiinctional when performing the variational calcu- 
lus. Actually, the conservation of the phase-space volume of fine- 
gra ined distr ibution is a natural consequence of the CBE. As proved 
by ( ICubarsill201Ql . the series of moment equations is equivalent to 
the CBE. So, if moment equations are included, the phase-space 
volume conservation is automatically satisfied. 

However, it would be very inconvenient to calculate the vari- 
ation using these moment equations directly. Like the example 
we will present in the next subsection with a truncated DF to the 
second-order expansion of velocity, we will use the equivalent se- 
ries of generalized virialization relations (equations |9], instead of 
the original moment equations dAl It . as the constraints of mechan- 
ical/virial equilibria. As a result, such a replacement will make the 
variational calculus more easy to perform. 

Such an approach to the statistical mechan ics for self- 
ravita ting systems stems from the work of Wh ite & NaravanI 
19871) . who applied the entropy of an ideal gas, = 
Anr^ pln(j>^^^ p~^^^), to gravitating systems, with the con- 
strained entropy, as in the following (their equation 12): 



St + pMt + \Et + MR) 



/dp 

\dr 



J 



(26) 



Here, fpiR) represents the truncation at the radius R. Notice 
that, besides the usual co nstraints of mass a nd en ergy conserva- 
tion (i.e. p.Mt and XEt). ' White & NaravanI lll987h also included 
the isotropic Jeans equation (i.e. the last integral term) as an addi- 
tional constraint. By extremizing Q with respect to p and p, they 
derived an equation for the mass density of fully relaxed galaxies. 

Although they did not obtain acceptable results, this approach 
still provides a coiTect line of thought, because, as we have men- 
tioned previously, the final relaxed states of self-gravitating systems 
are not only statistical equilibrium states, but also mechanical equi- 
librium states. Hence, the mechanical equilibrium equations should 
also be considered as constraints to the total entropy St, and the 
Jeans equation represents the second-order mode of the mechani- 
cal equilibri a of self-gravitating system s. Below, we generalize the 
approach of IWhite & NaravM (Il987h to formulate the methodol- 
ogy of the statistical mechanics for self-gravitating systems. 

With the DF F{x,v) given by equation l l22t . the total entropy 
from equation J20t is 



F{x,v) \nF{x,v)d^xd\ 



E 



Afe , 



J {x)viV^vsF{x, v)d^j:d^v 
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= ^2 J ^k,m,nix)Mk,m,nix)d'^X, (27) 

where we used equation ( |25t for the last equality. Because 
these moments, Mk,m.n{x), depend on the expansion coefficients, 
^k.m,nix), SO the total entropy Stot is a functional of \k.m,nix)- 

From equation i25\ . the mass density p{x) and kinetic energy 
density ek(j^) are 

p{x) = Mo,o,o(x), 

6kW = l-p(x)(vl(x)+^{x)+4{x)) 



{M2,Qfi{x) + Mo,2,q{x) + Mo,o,2(a:)). 



(28) 



So the total mass, kinetic, and potential energy of the system are, 
respectively. 



Mt 



E-K 



p{x)A X, 

e]^{x)d^x, 

Gp{x)p{x'] 
2\x-x'\ 



j3 j3 / 



(29) 



Here, the integrations are performed over all the spatial volume, 
and hence the total energy is Et = -Ek + Ey- As is the case for 
usual thermodynamic systems, in order to derive the proper ther- 
modynamic relationships of self-gravitating systems using the en- 
tropy principle, the two common constraints of mass and energy 
conservation should be applied. 

As addressed in Section 12.21 the series of moment equa- 
tions dAl lb does not constitute a closed equation system, because 
there are no equations for pv^, the one-dimensional transverse mo- 
ments of all orders {k — 0, 2, 4, . . .). This implies that Xk,m,n{x), 
the coefficients of the Taylor expansion in equation i22\ . cannot be 
completely determined (see Fig. [T}. It is the task of statistical me- 
chanics to determine these coefficients completely. In Table [T] we 
summarize the numbers up to the lOth-order for both the moment 
equations and the expansion coefficients. 

We have mentioned previously that mechanical equilibria 
must play crucial roles in establishing the statistical equilibrium 
states of self-gravitating systems. The moment equations dAI It . or 
their equivalent virialization forms (equations |9), represent the me- 
chanical/virial equilibria of various modes in spherically symmetric 
systems. We formally denote these virialization equations as non- 
spherical 'Eqns(x, T, W, P) = 0'. Therefore, the constrained en- 
tropy is 



+ 



E 

k.m.n 



PiMt - 

fk,n 



P2Et + 

,„{x)Eqns{x,T,W,P)d''x, (30) 



where pi, p2 and fk.m.n (x) are Lagrangian multip liers. Compared 
with the oris inal form of IWhite & NaravanI (Il987h . it can be seen 
that the mechanical/viral equilibria operate in very complicated 
ways, by including all orders of virialization equations into the 
constrained entropy. Because Stot, A/t and Et are also function- 
als of the velocity moments, by performing the variational calculus, 
5Stot,c ~ 0, with respect to these moments, Mk,m,n (x) [or equiva- 
lently \k,m,n (x)], we can formally obtain the differential equations 
concerning Mk,m,n{x) [or Xk,m,n{x)], which characterize the sta- 
tistical equilibria of the systems. These statistical equilibrium equa- 
tions, together with the mechanical equilibrium equations, provide 



Table 1. Numbers of moment equations and Taylor expansion coefficients 
corresponding to the order of moments (see equations lAll l and 1221 . Col- 
umn 1 is the order of the velocity moments, column 2 is the number of 
moment equations of the order, and column 3 lists the cumulative number 
of the moment equations up to this order. Column 4 is the number of Taylor 
expansion coefficients of the order and column 5 is the cumulative number 
of expansion coefficients up to the order. 



Order 


Velocity 


Cumulative 


Taylor 


Cumulative 


of 


moment 


moment 


expansion 


expansion 


moments 


equations 


equations 


coefficients 


coefficients 











1 


1 


2 


1 


1 


2 


3 


4 


2 


3 


3 


6 


6 


3 


6 


4 


10 


8 


4 


10 


5 


15 


10 


5 


15 


6 


21 



a complete description of the equilibrium states for self-gravitating 
systems. 

Finally, we emphasis that we should not directly include the 
steady-state CBE into equation i30\ as the constraint in place of the 
series of virialization equations, because the entropy Stot is con- 
structed with the coarse-grained DF F{x,v), which does not satisfy 
the CBE. Moreover, there is no simple relation between F{x, v) and 
the fine-grained DF f{x,v). 

We further explain the approach to the statistical mechanics 
that we have developed in the following example with a truncated 
DE 



4.3 Truncated DF of self-gravitating systems 

For the spherically symmetric systems, we truncate the DF F(x, v) 
in equation i22\ to the second order, and we directly replace the 
expansion coefficients by the corresponding moments as 



F{r, Vr,vo,v^) 



i2n) 



-exp(- 



2^ 



2a? 



2a?' 



(31) 



Here, the density p and the velocity dispersions a^ and a^ are all 
functions of r. Because of the spherical symmetry, the dispersions 
of Vff and are the same, so we uniformly use to represent 
Vg or v'^. As a result, the total (Boltzmann-Gibbs) entropy of the 
system with this truncated DF is 



Sbc 



Flni^d'^jcd-'v 



In 



p{r) 



ar{r)a?{r) 



In 27r) 



p{r)4:Tvr dr. 



(32) 



If we drop the unimportant constant term, (3/2) (1 + ln27r), and 
extract the specific entropy from equation (1321 . we have 



In 



= In 



1/2 
Pr Pt 

o5/2 



(33) 



where pr = pa,, and pt = pat are the radial and one- 
dimensional tangential pressure, respectively. It can be seen that 
the specific entropy of equation ( I33t greatly resembles that taken by 
White & Na rayan ( 1987), except for the anisotropic velocity distri- 
bution (i.e. Pr 7^ Pt). 

The Jeans equation l|3j can be transformed into the virializa- 
tion form of equation JI It 
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2fetot(r) + u(r) = Aur Pr{r), 

where ktot{r) and u{r) are the kinetic and potential energy con- 
tained in the r-sphere of the system. To be concise, fctot (r) is ex- 
pressed as 



fctot (r) = kr{r) + 2kt{r), 
with (see eguationfTOll 



(34) 



fcr(r-) = T2,0 = - 



kt{r) 



To,: 



47rr'^Pt(r')dr', 



(35) 



are the radial and one-dimensional tangential kinetic energy within 
the r-sphere, respectively. The potential energy u(r) is expressed 
as 



u[r) = Wo,o 



'sir 
-AtvG 



2\r' — r"\ 



m{r')p{r')r'dr' , 



(36) 



where the integral is confined within the volume of the r-sphere. 
Sir, and m(r) is the mass function 



m(r) 



47rr'^p(r')dr'. 



(37) 



With these variable definitions, the total mass A/t, the total 
kinetic energy Ek and the total potential energy E\r of the system 
are, respectively. 



Mt = m(r) 



Ek = fctot (r) 



Ev — u(r) 



.(38) 



From equations ( I35t and ( I37t . we perform the following vari- 
able transformation by calculating the first derivatives as 



kl 



From equation ( I36t . we have 

, , ^ Gmm' 
u — —AirGmpr = . 



Pt 



K ^ 



(39) 



(40) 



The prime in equations ([39} and ( I40t indicates the first derivative 
with respect to the radius r. 

Therefore, the total entropy of the system from equation i 

is 

fcr ^ (fctot fcr) 



St = 



+ ln(2^/Vr^) \ dr. (41) 



With the variable transformations of equations (I34t . ( I36t and BTb . 
the mass and energy conservation constraints for the variational cal- 
culus are converted in order to satisfy the fixed endpoint conditions 
of equation B8t . 

With equation ([39}, equation dllt can be re-expressed as 

2fctot(r) + u{r) - 2rfc;(r) = 0. (42) 

We use this relation, instead of the original Jeans equation ([3]l, as 
the constraint of the mechanical equilibrium. As a result, such a 
transfoiTnation will make the variational calculus more easy to per- 
form. Additionally, the mathematical relation between u and m, 
m' in equation \40h provides another additional constraint. 

With all these constraints considered, the total constrained en- 
tropy of the system is 



St. 



H{r,m, m, fctot, fctot, fcr, fcr, u, u)dr 



dr < m' In 



fcr ^ (fctot fcr) 



^/5/2 



I / 1 /r»5/2 

+ m ln(2 ' 7rr 



+ /i(r)(2fctot +u- 2rfc;) + f2{r){u' + ^^^) 



.(43) 



Here, H denotes the integrand (i.e. the terms enclosed in the braces) 
and /i and /2 are two Lagrangian multipliers. Then, performing 
the standard variational calculus, 5St,c = 0, with respect to the 
variable pairs (m, m'), (fctot, fctot), (fcr, fcr) and (u, u), we obtain 
the following equations: 



1 dlnpr dlnpt 

2 dr dr 

dr Pt 

= 2/ir + Ci, 



5 d In p , Gm 



Gm 



2 dr 



+ f2 = f2^^, 



Pr 



Pt 



d,f2 , 

d7^^^- 



(44) 
(45) 
(46) 

(47) 



Here, again, the prime in equation ( I44t indicates the first derivative 
with respect to r. In equation J46b . Ci is an integration constant, 
but the end point of fcr(r) at r — >■ oo is not fixed. So, from the 
variable end-point condition, we have 



dH 

dk'r 



P P Of 

2/ir 

Pr Pt 



= 0. 



From equation ( I46t , we can determine Ci 
grangian multiplier /i is 



/i = 



2rpr 



P 
2rpt 



(48) 

0. Thus, the La- 
(49) 



From equations ( 145 l l and ( I47t . we can obtain the other Lagrangian 
multiplier /2 as 



f - P 
2pt 



A, 



(50) 



where A is an integration constant. From equations ( 145) and ( 149) , 
we have 



^(P_) = J. fL 

dr Pt rpr rpt ' 
which can be further expressed as 
dlnp dlnp 
dr dr rpr 



r 



(51) 



(52) 



If we substitute the two Lagrangian multipliers into equation J44b . 
and use the Jeans equation ^ and equation i52\ . we can simplify 
equation ( 1441 ) as 



3 dlnp 
2 dr 



1 dpr 
Pt dr 



2pr 
rpt 



2 _ ^Gm 



(53) 



If the virialization relation, equation dllt or (142) . is not included 
as an extra constraint, then in this case /i — 0. So, /2 is just a 
constant, and from equation ([46} or equation ( 149) , we have pr — 
Pt- Hence, equation ( 1441 ) degenerates to 



3 dlnpr 
2 dr 



5 dlnp 
2 dr 



^A 



Gm 



(54) 



which is exactly th e equation that we derived in a previous study 
dHe & Kan j bold) . According to our latest results, this equation 
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should not be accurate. It can work well only in combination with 
the Jeans equation, and also with an empirical density profile (e.g. 
the Einasto form). 

To summarize, with the truncated DF of equation ( 131) , the 
complete set of variables contains only the density profile p{r) and 
the radial and one-dimensional tangential pressure pr (r) and pt{r). 
The Jeans equation ([3) and equations l |52t and l |53t , derived from 
our statistical-mechanical approach, constitute the complete equa- 
tion system for the three variables p, pr and pt. As analysed in 
Appendix [A] a separable DF such as equation ( |3U will not yield 
correct results. However, if the theoretical framework of the statis- 
tical mechanics that we have formulated is correct, we can expect 
that the predictions will be more and more accurate compared with 
the simulation results, as more and more expansion coefficients of 
F{x,v) in equation i22l are included. 

We see that, in our statistical-mechanical approach, as more 
and more moments and moment equations are involved, more 
and more expansion coefficients are needed to express the DF of 
equation i22\ . Thus, we can obtain any possible velocity distri- 
bution with the Boltzmann-Gibbs statistics. Hence, the idea that 
Boltzmann-Gibbs statistics can only produce a canonical distri- 
bution, s uch as a Maxw ell-Boltzmann distribution, is not correct 
(cf. H ansen et all 120051) . Kafri (20 09) presented a comment on 
non-extensive statistical mechanics ( iTsallisI 19881) . and pointed out 
that the Boltzmann-Gibbs statistics can also yield a long tail (i.e. 
power law) distribution. Additionally, it might also be unneces- 
sary to invoke the general - function o ther than the Boltzmann- 
Gi bbs entropy (cf. Tremaine et al. See also the discussion 

bv lBouchet & BarrdI ( I2OO6I) about the validity of Boltzmann-Gibbs 
statistics when using large deviation tools. 



5 DISCUSSIONS AND CONCLUSIONS 

The statistical mechanics of self-gravitating systems is a long- 
standing puzzle, which has not yet been successfully solved. In this 
paper, we have formulated a systematic theoretical framework of 
the statistical mechanics of spherically symmetric self-gravitating 
systems, using approaches that differ significantly from the conven- 
tional statistical mechanics of short-range interaction systems. 

The equilibrium states of self-gravitating systems also consist 
of mechanical equilibria. By multiplying the steady-state CBE with 
various powers of velocity, and integrating over all velocities, we 
obtain the moment equations of various orders. In Appendix lAl we 
demonstrate that the DF / is inseparable. As a result, the moments 
are infinitely many and independent of each other. These moments 
constitute the endless moment equations dAllI ). which character- 
ize the great complexity of the mechanical equilibria. The moment 
equations can be further transformed into the generalized virial- 
ization equations (|9), and hence mechanical equilibria are actually 
identical to the generalized virial equilibria. 

One subtlety that should be pointed out is that the DF /, which 
satisfies the CBE, is the fine-grained DF. All the velocity moments 
are derived from this fine-grained DF, whereas the important DF 
for statistical mechanics is the coarse-grained DF. Although the 
coarse-grained DF does not obey the CBE, fortunately, we find that 
the coarse-grained moments differ from the fine-grained moments 
only by up to second-order small quantities determined by the size 
of the finite volume element. So, we can safely use these moment 
equations to characterize the mechanical equilibria of the system. 

Ergodicity, or the equiprobability principle, is fundamental for 
equilibrium statistical mechanics. However, ergodicity is invalid for 



the whole coUisionless self-gravitating system, because, for exam- 
ple, a gravitating particle moves in a static gravitational field along 
a specific curve, on which the energy and angular momentum of 
the particle are conserved. Such a regular motion of the particle is 
definitely not an ergodic motion, in that the particle cannot traverse 
any place in the system, and also cannot attain any value of veloc- 
ity. Even if the particles are treated as indistinguishable, the validity 
of global ergodicity cannot be recovered, either. According to the 
ensemble theory, the probability distribution of the microstates of a 
closed system at a fixed energy is derived by assuming that all states 
on the energy hypersurface in phase space have equal probability. 
However, the hypersurface of a constant energy for an isolated self- 
gravitating system is unbounded. Thus, the microcanonical proba- 
bility distribution cannot be defined, which implies that ergodicity, 
or equiprobability, is broken. 

It seems to be a disaster if global ergodicity is invalid, because 
the equilibrium statistical mechanics cannot exist without the va- 
lidity of ergodicity. Fortunately, ergodicity can be re-established 
locally in the following manner. Within an imagined finitely small 
volume element of the system, none of the coUisionless gravitat- 
ing particles can be tracked, because these particles are not con- 
fined within this limited volume element. Hence, these particles 
should be treated as indistinguishable. At some instant, these par- 
ticles, randomly and homogeneously populated in the element, ex- 
hibit a static distribution pattern. These different patterns of differ- 
ent instants naturally constitute an ensemble. When the system is 
in its equilibrium states, the background field quantities (i.e. all the 
velocity moments) are local invariants, which define a hypersur- 
face for this volume element. It is reasonable to assume that every 
member of the ensemble of the local distribution patterns occurs 
with equal probability in this hypersurface. Intuitively, if within 
the finitely small element, the particles are randomly and homoge- 
nously distributed, and the velocity distribution is independent of 
spatial coordinates, then the local ergodicity should be valid and 
can be re-established. 

We assume that the minimum phase-space element (i.e. the 
microcell) can accommodate only one gravitating particle, as it is 
impossible for two gravitating particles to be close to each other 
but simultaneously to be relatively static. Because of the local er- 
godicity, the particles will populate these microcells with equal 
probability, which also indicates that all the microstates for this 
element occur with equal probability. So, the microstates of the 
indistinguishable gravitating particles within this volume element 
should be counted with Fermi-Dirac-like statistics. Moreover, be- 
cause both the direct gravitational interaction and the statistical 
correlation between any two gravitating particles can be neglected, 
which is exactly the same reason to justify the CBE, it is reason- 
able to assert that the microstates of one volume element do not 
depend on those of any other elements (i.e. local microstates are 
spatially independent). Hence, the total microstate number of the 
whole system can be expressed as the product of the microstate 
number of all the elements. Based on the local ergodicity, and by 
using Fermi-Dirac-like statistics, with a non-degenerate condition 
and the spatial independence of the local microstates, we rederive 
the Boltzmann-Gibbs entropy, which is exactly the correct entropy 
form suitable for coUisionless self-gravitating systems. Finally, the 
Boltzmann-Gibbs entropy is additive and h ence extensive, which 
means that the non-extensive Tsallis' ('1988*) entropy as well as the 
general //-function (Tremaine et al. 1986) are not necessary. 

The equilibrium states of self-gravitating systems consist of 
both mechanical/virial and statistical equilibria, of which the lat- 
ter should be characterized by the statistical equilibrium equa- 
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tions, derived from entropy principle. Apart from the usual con- 
straints of mass and energy conservation, we demonstrate that the 
series of moment equations dAl it or their equivalent virialization 
equations ^ are not only employed to characterize the mechan- 
ical/virial equilibria, but also act as additional constraints on the 
system's entropy functional when performing the variational calcu- 
lus. This approach is an extension of the original prescription by 
IWhite & NM-avanl ( Il987h . who just used the isotropic Jeans equa- 
tion as an extra constraint. 

We further illustrate this statistical mechanical approach with 
a truncated, separable DF. In this case, the complete set of vari- 
ables can be chosen to be the density profile and the radial and 
one-dimensional tangential pressures. The two equations (152b and 
( I53t , derived from the entropy principle, together with the Jeans 
equation l[3), constitute the complete equation system for the three 
variables. As analyzed in Appendix|A] such a separable DF will not 
yield correct results. However, if the theoretical framework of the 
statistical mechanics that we have formulated is correct, we can ex- 
pect the predictions to be more and more accurate, compared with 
simulation results, because more and more expansion coefficients 
of the DF are included. As a result, we can produce any possible 
velocity distribution using the Boltzmann-Gibbs statistics. Hence, 
the idea that Boltzmann-Gibbs statistics can only produce a canon- 
ical distribution, such as a Maxwell-Boltzmann distribution, is not 
correct. 

We also discuss the problems of negative specific heat and en- 
semble inequivalence. We propose that these problems might origi- 
nate from the misuse of ensemble theory. This is because ensemble 
theory is established only upon the validity of global ergodicity, 
which is not valid for self-gravitating systems. As a result, the re- 
lationships about the heat capacity of equations ( |13t and l ll4t are 
rootless for self-gravitating systems. Without the validity of ensem- 
ble theory, these two relationships should be, at most, considered as 
two different definitions of the 'heat capacity' for self-gravitating 
systems. 

We have formulated a systematic theoretical framework of 
the statistical mechanics of spherically symmetric self-gravitating 
systems, within which we encounter many difficult mathematical 
problems, such as the evaluation of both the entropy functional and 
the velocity moments, and the variational calculus of deriving the 
equations of statistical equilibria. These will be the tasks of future 
investigations. 
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APPENDIX A: COLLISIONLESS BOLTZMANN 
EQUATION, VELOCITY MOMENTS AND MOMENT 
EQUATIONS 

The CBE in spherical coordinates (r, 6, (j)) is jBinnev & Tremain3 
12008.) 
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df .df -df -df . df . df . df „ 

at or do ocp ovr ovg ov^ 

where the time derivatives of the coordinates can be expressed in 
terms of the velocity components: 

r r sm o 



(A2) 



The Lagrange equations give the components of the acceleration: 

Vg + vl a$ 



Vr ~ 



r dr ' 

v'i cot 6 — VrVe 1 9$ 



V0 = 



Vtf, = 



r r 86 ' 

-v^Vr — v^vg cot 6 1 9$ 



r sm t 



(A3) 



For a spherically symmetric system in its steady state, all the terms 
in equation JAU . including d/dt, 8/89 and 8/8(j), should be van- 
ishing. Furthermore, 8^/89 and 8^/8<ji in equation l lA3b should 
also be zero. 

With all these to hand, we can multiply the steady-state CBE 
by various powers of velocity, and integrate over all velocities, to 
obtain a series of moment equations. 

For a spherical system, the fine-grained DF /(x, v) reduces 
to the form of /(r, n,., ve,V4,), where Vr is the radial velocity and 
ve and are two tangential velocities. The various orders of ve- 
locity moments, such as p(r), pv", pvg, pv^, pu"Ug", pu" «™ and 
pv'gV'T, can be uniformly defined as 



pVrV^'V'^ 



(A4) 



Here, k, m, n are non-negative even integers (i.e. k,m,n — 
0,2,4,...). 

If we multiply the steady-state CBE by v'^ ^w^^S and inte- 



grate over all velocities, we obtain 

(m + n + 2) 



pVrV^V^ 



pVr Vg ^ V" 



pVr Vg V , 



mcot9 ™_i „+2 , (w + 1) cot 9 „+i 

PVr Vg V^^ + pVr Vg ^ V^ 

N d$ — j-^ 

= -ik-i)^pv^-'v^v;, 

which can be further divided into two relations. The first relation is 

(m + n + 2) 



(A5) 



— ipv^v^v^) + 



pv^v^v^ 



-pVr Vg V'^ 



{k - 1) -J-T 



-pVr VgV^ 



(n+ l)(n- 1) „,+3 2 _ 

m(m + 2) " 

(m-2)!!(n + l)!! 

(m + n)!! « 



(A8) 



(n + 1) 



fc-l m-l ,n+2 
Vr Vr, V, 



Am - 2) -J— I 



(n+ l)(n + 3) 
m!!(n- 1)!!- 



Ug — 



(m + n)!! ^ 



(A9) 



Equations dASI l and l |A9l > indicate that 



fc — 1 m — 1 n + 2 



«fl V 



y^^w^^^'U^ are not independent. From equations (IA7) . ( IA8t and 
( |A9t , with the replacement fc — 1 — > fc, and m + n + 1 — m, we 
can obtain 



(AlO) 



where fc, m are even numbers (fc, m — 0,2,4,.. .). With equa- 
tions (|A7} - dAlOb , and with the replacement fc — !> fc + 2, and 
m + n m, equation iA6\ can be reduced to 



(m + 2) —2 



-^{pvr vY") + ■ ^ ■ pVr ■ -vi- 



(fc + 1) (m + 2) ^ ,.fe^„,+2 



d$ 



r (m + 1) 



pVrV^ 



^ — ik + l)——pVrV 

dr 



(All) 



where is either vg or w^, with fc,m = 0,2,4,.... Gener- 
ally, there are A'^ independent 2A^th-order moment equations, with 
different combinations of fc a nd m in e quation dAl It . See also 
iDeionghe & Merritll ( Il992h and|A^ ( l201ll) for the derivation of this 
equation. 

We would like to investigate whether these moment equations 
can be reduced to their counterparts of lower orders. If the DF is 
uncorrelative, i.e. / can be separated as 



f{Vr,V0,V^) = h(Vr)f2{vg)h{v^), 



(A 12) 



where /i, /2 and fz are one-dimensional DFs, then the velocity 
moments of equation dA4b can also be separated as 



Ur Ug u^. 



Setting n = 0, and with the replacement m — )■ 2rn + 1, from 
equation dA7t we obtain 



(2m + vl 



= {2m + l){2m^iyg--'{vl) 



2\2 



, d<E> 



= -(k-i)^pvr'v^v;, 



(A6) 



where fc should be an even number. In this case, n and m must also 
be even numbers, with fc ^ 2, and m, n ^ 0. The second relation 
is the following recursive relation 



mvr\^-\2 = (n+ l)wr 



fe-l„m+l. 



(A7) 



where fc, m are odd numbers, with fc, m 1, and n is an even 
number, with n ^ 0. Using equation dA7t repeatedly, we have 



= (2m + l)!!(«|)™+^ 
= (2m + l)!!(^r+\ 
Generally, with tit being either vg or w^, we have 



^2m+2 ^ (2m+ l)!!(^;?)™+\ 



(A13) 



(A14) 



which means that the moments of Vt are not independent. In this 
case, using the Jeans equation the 4th-order moment equa- 
tion ^ can be reduced to 



^2 d ( 2)+ 2 
dr r 



(A15) 
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This can be re-expressed as 




(A16) 



where /3 is the ani sotro py parameter. From the simulations by 



creases with the radius r from zero to the maximum, and then de- 
creases to about zero again at the outskirts of the dark halos (i.e. 
/9 ^ for most of the radii, except for the possible negative value 
at the outer boundary). Hence, nf resulting from equation ( IA16) 
should monotonically decrease with r for m ost of the radii, which 
is in contrast to the same simulation results of NavaiTO et al. ( 2010) 
and iLudlow et alj ( I2OIO,) . This suggests that equation 1 IAI6I 1 and 
equations ( IA13t - dAlSb are not correct. Hence, the DF / is in- 
separable, which implies that these moment equations cannot be 
reduced to the lower-order counterparts. 
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